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The concept of heavy-particle radioactivity (HPR) is changed to allow emitted particles with 
Ze > 28 from parents with Z > 110 and daughter around ^"^Pb. Calculations for superheavy (SH) 
nuclei with Z=104-124 are showing a trend toward shorter half-lives and larger branching ratio 
relative to a decay for heavier SHs. It is possible to find regions in which HPR is stronger than 
alpha decay. The new mass table AMEll and the theoretical KTUY05 and FRDM95 masses are 
used to determine the released energy. For 124 we found isotopes with half-lives in the range of ns 
to ps. 
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In recent years the heaviest elements with atomic num- 
bers up to Z = 118 have been synthcsised [l[ either 
with cold fusion reactions having the ^°®Pb or ^"^Bi tar- 
get [1, Q or with hot fusion induced by ^^Ca projectiles 
[J,Q. Attempts to produce Z = 120 are reported [y] and 
new experiments are presently running at GSI Darmstadt 
0. The main experimental difficulty in identifying the 
new superheavy (SH) elements is the low probability of 
their formation, and the separation of the short lived 
compound nucleus from the very high flux of incident 
projectile nuclei. The lowest cross-section of 55 fb was 
measured at RIKEN Q where one decay chain of ^^^113 
was observed during 79 days with beam of ^"Zn on ^"^Bi 
target. After naming copernicium, Cn, Z ~ 112 sug- 
gested by GSI scientists, lUPAC recommends that the 
Dubna-Livermore collaboration be credited with discov- 
ery of new elements 114 and 116. 

It is generally agreed that the term SH element, intro- 
duced [8| in 1958, is a synonym for elements which exist 
solely due to their nuclear shell effects. The lightest SH 
is Z = 104 Rf with half-lives of different isotopes around 
1 min. This is 16 orders of magnitude longer than the 
expected nuclear lifetime of 10"^'^ s these isotopes would 
survive without any shell stabilisation. Spontaneous fis- 
sion, the dominating decay mode in the region around 
Rf, becomes a relatively weaker branch compared to a- 
decay for the majority of recently discovered proton-rich 
nuclides. Extensive calculations of fission barriers and 
half- lives have been published @. 

Despite the important experimental and theoretical de- 
velopment there are still several unanswered questions 
related to the magic numbers, production cross sections, 
and decay modes. Besides beta decay, only alpha decay 
and spontaneous fission of SH nuclei have been exper- 
imentally observed up to now. We would like to take 
also into account heavy-particle radioactivities (HPR) 



Since 1984 [l2|, the following HPR have been exper- 



imentally confirmed 13| in heavy parent nuclei with 
Z = 87~ 96: "c, 209, 23p^ 22,24-26^6, 28.30^5, 32.343^ 

with half-lives in good agreement with predicted values 
within analytical superasymmetric fission (ASAF) model 
(see the review [TJ] and references therein). Almost al- 
ways the corresponding daughter nucleus was the doubly 
magic |2^Pbi26 or one of its neighbors. The newest mea- 
surement of ^^C decay of 223 Ac [l3| was one of the pos- 
sible candidates for future experiments mentioned in the 
systematics [3l showing that the strong shell effect due 
to magic number of neutrons, = 126, and protons, 
Zd = 82, present in order to lead to shorter half- lives 
was not entirely exploited. 

The shortest half-life of Tc = 10^^ °^ s corresponds to 
^^C radioactivity of ^^^Ka and the largest branching ratio 
relative to alpha decay ba = Ta/Tc = 10~®'^ was mea- 
sured for the ^■^C radioactivity of ^^^Ka,. Consequently 
HPR in the region of transfrancium nuclei is a rare phe- 
nomenon in a strong background of a particles. Several 
attempts to detect ^^C radioactivity of the neutrondefi- 
cient ^^*Ba with a daughter in the neighbourhood of the 
double magic gg'^Snso, predicted to have a larger ba, have 
failed. 

In order to check the possibility of extrapolations from 
Ae = 14 — 34 emitted clusters already measured in the 
region of emitters with Z — 87 — 96 to SHs up to 
124, where one may find an emitted particle as heavy 
as ^^^Mo, we estimated within ASAF model the half- 
life for i28Sn emission from 256Fm (Q = 252.129 MeV) 
and for ^^OTe emission from ^^^Rf (Q = 274.926 MeV): 
log^oTpmis) = 4.88 and logioTfl/(s) = 0.53. They are 
in agreement with experimental values for spontaneous 
fission [13: 4.02 and 0.32, respectively. 

There are many other theoretical approaches of the 
HPR e.g. Refs. 3"21|. Any spontaneous emission of a 
charged particle from atomic nucleus may be explained as 
a quantum mechanical tunnelling of a preformed cluster 
at the nuclear surface through the potential barrier [2^ . 



2 



Microscopic calculations of cluster formation probability 
and of barrier penetrability have been performed 2^, 21 1 
by using the R-matrix description of the process. The 
half-life, Tc, is expressed as 



h\n2 



(1) 



where Tc is the decay width and h is the Plank constant. 
A universal decay law for a emission and HPR was re- 
cently developed [2l| based on this theory. 

We should change the concept of HPR, previously 
associated to a maximum Z^"'^ = 28, allowing to pre- 
serve its main characteristics in the regions of SH with 
Z > 110 i.e. in a systematic search for HPR we shall con- 
sider not only the emitted particles with atomic numbers 
2 < Ze < 29, as in previous calculations, but also heavier 
ones up to Z™"^ = Z — 82, allowing to get for Z > 110 an 
atomic number of the most probable emitted HP Z^ > 28 
and a doubly magic daughter around ^°*Pb. 

Calculations are performed within the ASAF model, 
very useful for the high number of combinations parent- 
emitted cluster (more than 10^) in order to check the 
metastability of more than 2000 parent nuclides with 
measured masses against many possible decay modes. 
We started with Myers- Swiatecki liquid drop model 
(LDM) [2^ adjusted with a phenomenological correction 
accounting for the known overestimation of the barrier 
height and for the shell and pairing effects in the spirit 
of Strutinsky method. 

The half-life of a parent nucleus AZ against the split 
into a HP or an emitted cluster A^Zg and a daughter 
AdZd is given by 



T = [{hln2)/{2E„)]exp{Koy + K^) 



(2) 



and is calculated by using the Wentzel-Kramers- 
Brillouin (WKB) quasiclassical approximation, according 
to which the action integral is expressed as 



/i = - / yj2B{R)E{R)dR 



(3) 



with B = fj. the reduced mass, K = Kqv + Ks, and the 
E{R) potential energy replaced by [E{R) — Ecorr] — Q- 
Ecorr is a correction energy similar to the Strutinsky [25| 
shell correction, also taking into account the fact that 
LDM overestimates fission barrier heights, and the ef- 
fective inertia in the overlapping region is different from 
the reduced mass. Ra and Rt are the turning points 
of the WKB integral. The two terms of the action in- 
tegral K, corresponding to the overlapping (Kov) and 
separated (Kg) fragments, arc calculated by analytical 



formulas [14 1 



Half-life calculations arc very sensitive to the released 
energy (Q value) obtained as a difference of the parent 
and the two decay product masses 



in units of energy. Even with the newly released ta- 
bles of experimental masses, atomic mass evaluation 2011 
(AMEll) [H] as a preview for the AME13 pubhcation, 
many masses are still not available for new SH, hence we 
shall use not only these updated tables for 3290 nuclides 
(2377 measured and 913 from the systematics) ending up 
a.t Z = 118 but also some calculated masses, e.g., Koura- 
Tachibana-Uno-Yamada (KTUY05) ^ and the finite- 
range droplet model (FRDM95) [H] with 9441 and 8979 
masses, respectively. 

In a systematic search for HPR we calculate with the 
ASAF model for every parent nucleus "^Z the half- lives 
of all combinations of pairs of fragments "^^Ze, ^''Zd 
with 2 < Ze < z™'^^ conserving the hadron numbers 
Ze + Zd = Z and A^ + Ad = A. Let us start with the 
results obtained by using the AMEll mass tables. An 
example of the time spectra obtained for different clus- 
ters emitted from the parent nuclei ^^^Ra and ^^^114 is 
shown in Fig. [ijversus the mass numbers of the light frag- 
ment. The symbols of the emitted HPR are given on the 
figure's legend. 
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FIG. 1. (Color online) Time spectra of different cluster emis- 
sions from ^^^Ra (left panel) and from the superheavy nucleus 
^**114 (right panel). The most probable emitted clusters from 
^^^Ra and 2*** 114 are and '^'^Ge, respectively, both leading 



to 



208 



Pb daughter nucleus. 



M - (A4 + Md) 



(4) 



From the left panel of this figure one can see that the 
shortest half-lives of ^^^Ra correspond to a decay and 
I^Cs radioactivity, respectively. Both these decay modes 
have been experimentally observed and there is a good 
agreement between the calculated values and measured 
data. Other HPR with half-lives < 10^^ s are: ^Be; 
I2,i3(^. i5]\j g^jjjj 16q — much longer half-lives. 

Similarly, on the right hand side of the figure we show 
calculated results for the SH nucleus ^^^114. Again a de- 
cay is the strongest decay mode and there is a good agree- 
ment between our calculations and the experimentally 
observed half-life T^. The time spectrum in the region 
of mass numbers of emitted particles around A^ = 80 is 
more complex looking similar to a fission fragment spec- 
trum. There are many HPR with Ze ~ 28 ~ 38 having 
Tc < 10^^ s. For the sake of clarity we only plotted the 
results corresponding to even-even emitted HP which are 
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FIG. 2. (Color online) Chart of heavy and superheavy cluster emitters with atomic numbers Z = 94 — 118. The Q values are 
calculated using the AMEll mass tables [2^. Black squares mark the Green approximation of the line of beta stability. One 
most probable emitted cluster is given for every parent nucleus. 



leading to shorter half-lives in the same way the ^^C ra- 
dioactivity of ^^^Ra is less probable than both ^^C and 
^^C spontaneous emissions. The most probable emitted 
HP from ^^^114 is 32Ge48 with a calculated branching 
ratio ba = 10^^ *'^. One should also take into account a 
competition of USeso with a magic number of neutrons 
Ne = 50 and a branching ratio ba = 10~^'^^. 

We proceed in a similar way with all parent nuclei with 
Z = 94 - 118 present on the AMEll mass table. The 
chart of cluster emitters from Fig. [2] is obtained by as- 
sociating to each parent only the most probable emitted 
cluster. The black squares mark the Green approxima- 
tion of the line of beta stability. All superheavy nuclei 
present on the AMEll mass table are proton-rich nu- 
clides with neutron numbers smaller than Njj on the line 
of beta stability. The experimentally determined ^I^S 
radioactivity of II^Pu, f^Si radioactivity of gt^Pu, and 
fISi radioactivity of g|^Cm are fairly well reproduced. 

New many types of HPR with > 28 may be seen 
on this chart: Cu, Zn, Ga, Ge, As, Se, Br and Kr. 
We used only one color for a given despite the fact 
that as the result of the calculations we obtained sev- 
eral isotopes of these elements, e.g. = 26, 28 for Mg; 
30, 32, 33, 34 for Si; 36, 38, 40, 41, 42 for S; 44, 46, 47, 48 for 
Ar; 48, 49, 50, 51, 52 for Ca; 50, 51, 52 for Sc; 53, 54, 55, 56 
for Ti; 57,58,59,60,61 for Cr; 60,62,63,64,66 for Fe; 
66, 68, 69, 70, 71, 72, 73 for Ni; 69, 71, 72, 73, 74, 75 for Cu; 
72, 74, 76, 77 for Zn; 75, 77, 78, 79 for Ga; 78, 80, 81 for Ge, 
81,83 for As; 82,84,85 for Se; 85,86 for Br, and 86,87 
for Kr. Only one mass value was obtained for the most 
probable emitted particles Be, C, Ne, P, CI, V, and Mn. 

As we previously observed [i^, many of the proton- 
rich SH nuclides are ^Be emitters, but they have a very 
low branching ratio ba- The general trend of a shorter 
half-life and a larger branching ratio when the atomic 
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FIG. 3. (Color online) Decimal logarithm of the half-lives 
of superheavy nuclei against cluster radioactivities versus the 
neutron number of the parent nucleus. Q values are calculated 
using the AMEll experimental mass tables [2^ (left panel) 
and the KTUY05 [13 calculations. 

and mass numbers of the parent nucleus increases may be 
seen on the left hand side of the figures |3] and SI obtained 
within ASAF model by using the AMEll mass tables to 
calculate the Q values. 

One can advance toward neutron-rich nuclei by us- 
ing the KTUY05 calculated mass tables, as shown in 
the right panels of these figures. When using KTUY05 
and FRDM95 masses for parent and daughter nuclei we 
take into account the nuclides stable against one proton, 
two protons, one neutron and two neutrons spontaneous 
emissions. If the calculated masses are reliable, then half- 
lives Tc in the range of nanoseconds to picoseconds for SH 
nuclei with Z = 124 (see the right hand side of Fig. [3| 
would make difficult or even impossible any identifica- 
tion measurement. More interesting for future experi- 
ments could be some even-even proton-rich isotopes of 
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FIG. 4. (Color online) Decimal logarithm of the branching 
ratio relative to a decay for cluster emission from superheavy 
nuclei versus the neutron number of the parent nucleus. Ver- 
tical dashed lines correspond to A'^ = 174, 186, 198. 



the 122 element with iV = 188 - 194 for which the neu- 
tron number of the Green approximation of the line of 
beta stability is Njs = 202. 

The pronounced minimum of the branching ratio at 
N = 186 in Fig.|3]is the result of the strong shell effect of 
the assumed magic number of neutrons N = 184 present 
in the KTUY05 masses. The half-life of a decay of a SH 
nucleus with A^ = 186 neutron number leading to a more 
stable daughter with magic neutron number Nd = 184 
is shorter by some orders of magnitude compared to the 
a decay of a SH with N = 184. Calculated branching 
ratios 6q, > 1 for Rf {Z = 104) only occur in very neutron- 
rich nuclei with N = 194 — 200 compared to = 166. 
Also their Tc half-life is extremely long. Similar results 
were obtained using the FRDM95 masses. 

In conclusion, the concept of HPR should be changed 
to allow spontaneous emission of heavy particles with 
atomic number larger than 28 from SHs with Z > 110 
and consequently daughter nuclei around the doubly 
magic ^°^Pb. The calculated half-lives Tc against HPR 
and the branching ratios relative to a decay ba are show- 
ing a trend toward shorter Tc and larger ba for heavier 
SH nuclei which are not synthesised until now. If the 
KTUY05 and FRDM95 masses used to calculate the re- 
leased energy Q are reliable, we expect to find for the 
element 124 many isotopes with half-lives in the range of 
nanoseconds to picoseconds, making practically impossi- 
ble to perform any identification experiment. Neverthe- 
less, there would be a chance to observe some proton-rich 
isotopes of 122 with branching ratios ba > 1- 

We are looking forward to receive experimental infor- 
mation about the decay modes of SHs with Z > 120, 
hoping to confirm the present calculations. There is also 
a need for developing more refined decay models as well 
as new calculated mass tables and new mass measure- 
ments. 
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